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Abstract 

We study the Ashtekar-Jacobson-Smolin equations that charac- 
terise four dimensional complex metrics with self-dual Riemann tensor. 
We find that we can characterise any self-dual metric by a function 
that satisfies a non-linear evolution equation, to which the general so- 
lution can be found iteratively. This formal solution depends on two 
arbitrary functions of three coordinates. We construct explicitly some 
families of solutions that depend on two free functions of two coordi- 
nates, included in which are the multi-centre metrics of Gibbons and 
Hawking. 
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1 Introduction 



In four dimensions the Hodge duality operation takes two forms to two forms. 
Given a four dimensional metric, the most important two form associated 
with it is the curvature two form R a b . It is therefore natural to be interested in 
four dimensional metrics whose curvature form obeys the self-duality relation 

R\ = *R a b , (1) 

where * is the Hodge duality operator. We will refer to such metrics as "self- 
dual" . Such metrics automatically have vanishing Ricci tensor, and so satisfy 
the vacuum Einstein equations with vanishing cosmological constant. Unfor- 
tunately, the only real Lorentzian self-dual metric is flat Minkowski space, 
so we choose instead to work with metrics with four complex dimensions. 

Physically these metrics may be of interest in attempts to quantise grav- 
ity, since they correspond to saddle points of the Einstein-Hilbert action, 
therefore giving large contributions to a path integral over euclidean met- 
rics Alternatively, it may be possible to interpret them as "one particle 
states" in a quantised gravity theory f2|. 

From a purely mathematical point of view these metrics are interesting 
since they are "hyperkahler". Hyperkahler manifolds are 4n dimensional 
manifolds (for n a positive integer) that admit a non-singular euclidean met- 
ric, g, with respect to which there exist three automorphisms, J 1 , of the 
tangent bundle which obey the quaternion algebra ||. In other words 

VJ* = 0, J*J J ' = + e^J*, (2) 

where V is the covariant derivative with respect to the metric g. In four 
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dimensions, it turns out that for a metric, g, to be hyperkahler it must have 
either self-dual, or anti-self-dual, curvature tensor ||. 

The problem of constructing metrics with self-dual curvature tensor has 
been tackled in several ways. The most direct approach is to formulate the 
problem in terms of partial differential equations || || . A more constructive 
approach is Penrose's 'Non-Linear Graviton' technique |§. Here, the task 
of solving partial differential equations is replaced by that of constructing 
deformed twistor spaces, and holomorphic lines on them. In practise this 
turns out to be just as difficult as solving partial differential equations, but 
in principle one can construct the general self-dual metric in this way. 

Here we concentrate on partial differential equations. We find a formu- 
lation which is similar to Plebanski's First Heavenly equation ||, but which 
can be viewed as simply an evolution equation. This means that the free 
functions in our solution are just a field and its time derivative on some ini- 
tial hypersurface i.e. two free functions of three coordinates. We construct, 
in a somewhat formal manner, the general solution to this equation. We also 
construct explicitly some infinite dimensional families of solutions to these 
equations. In the appendices, we show how this formulation is equivalent to 
Plebanski's. 

2 Construction of Self-Duality Condition 

In the equations for complex self-dual metrics were reformulated in terms 
of the new Hamiltonian variables for General Relativity introduced in ||. 
By fixing the four manifold to be of the form M. = Sxi? and using the 
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coordinate T to foliate the manifold, they reduced the problem of finding 
self-dual metrics to that of finding a triad of complex vectors {Vj : i = 1, 2, 3} 
that satisfy the equations 

DivVi = 0, (3) 
^V 4 = i e ^[V,,V fc ]. (4) 
Defining the densitised inverse three metric 

q ab = V«V J b 5 tJ , (5) 

we recover the undensitised inverse three metric q ab by the relation q ab = qq ab , 
where q = det q a b = (detq ah )~ 1 . If we now define the lapse function iV by 
iV = (det qab) 1 ^ 2 then we find that the metric defined by the line element 

ds 2 = N 2 dT 2 + q ab dx a dx b (6) 

is self-dual. 

Later, it was found that this triad of vectors could be related to the 
complex structures J* that hyperkahler metrics admit 0. Given a self-dual 
metric, we choose local coordinates (T, x a ) to put the line element in the 
form of equation (^). If we define the triad of vectors Vj = — J*(*, Dt), then 
these vectors will satisfy @ and (|J). 

Here we will concentrate on the problem of finding local solutions to equa- 
tions (Q) and (Q). We thus introduce a local coordinate chart (X,Y,Z) on 
the three surface, E, with its natural flat metric and connection. Thus (0) 
becomes just 

i) 



The crucial step is to realise that we can write equation (^) as 

[^,n = \^k[y 3 ,v k ]. (8) 

If we consider only euclidean metrics, then we take the V, to be real. In this 
case we define two complex vectors A, B by 

B = V 2 -*V 3 . (10) 
which, by virtue of obey the Lie bracket algebra 

[A, B] = 0, [A, B] = 0, [A, A] + [B, B] = 0, (11) 

where ~ denotes complex congugate. We can generalise these equations by 
considering four complex vectors U, V, W and X that satisfy the relations 

[U,V] = 0, (12) 

[W,X]=0, (13) 

[U,W] + [V,X] = 0. (14) 

Here we are thinking of W and X as "generalised complex conjugates" of 
U and V respectively. By Frobenius' theorem, we can use fll2"|) to define a 
set of coordinates (t, x) on the 2 (complex) dimensional surface defined by 
vectors U and V, and take U and V to be 

u = l' v = |- (15) 

We can now foliate our whole space using the coordinates (t,x,y,z). The 



equation (II) then becomes c^W + c^X = 0. This means there exists a 



vector field Y such that W = d x Y, X = —dtY. Thus we are only left with 
the problem of solving for vectors Y that satisfy [dtY, d x Y] = 00. We expand 
W and X as 

W = d t + f x d y + g x d z , (16) 

X = -f t d y -g t d z . (17) 

(The reason for the d t term in W is, as alluded to above, we are thinking 
of W as a sort of complex conjugate of U = d t . Although this argument 
only seems sensible for t a real coordinate, we are still perfectly at liberty to 
expand W in this way if t is complex.) If, by analogy with (0), we impose 
-£^W a = -r^X a = 0, then we find that there exists a function h(t,x,y,z) 



such that f = h z ,g = —h y . Imposing ([Of), we find that there exists a 
function a(t,x) such that 

ht + h xz h ty - h xy h tz = a(t,x). (18) 

We can absorb the arbitrary function a into the function h, and conclude 
that we can form a self-dual metric for any function h that satisfies 

h tt + Kzh t y - h xy h tz = 0. (19) 

This is just an evolution equation. Thus we can arbitrarily specify data h 
and h t on a t = constant hypersurface and propagate it throughout the space 
according to ( |T9"D to get a solution. For example, if we expand h around the 
2 It was only after this work was completed that I learned of jlO) where the ideas 
developed so far were found independently. From here onwards, however, our treatments 
are different. 
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t = hypersurface, and insist that h is regular on this surface, then h is of 
the form 

t 2 t 3 

h = a (x, y, z) + a x (x, y,z)t + a 2 (x, V,z)— + a 3 (x,y,z)— + ... (20) 

Substituting this into ( |19"D shows that olq and a% are arbitrary functions of 
x, y and z. a 2 , a 3 . . . are then completely determined for chosen a and a x 
by 

a 2 — dOxydlz — O'OxzO'ly^ (21) 
= O-Oxy d2z — dOxz^y + a lxy 0-1 z ~ °<lxz a ly, (22) 

and so on. Thus, in principle, we have a solution that depends on two 
arbitrary functions of three coordinates. It is interesting to compare our 
equation (|l^) with Plebanski's First Heavenly equation 

Qpq Qpq ^Ipp^lqq 1. (^^) 

Here it is not so obvious what our free functions are, and an expansion 
along the lines of ( |2"0"D doesn't work. (It is shown in the appendix how to 
get equation (^3j) from our equation, showing that the two approachs are 
equivalent. Thus for any self-dual metric there will exist a corresponding 
function h that satisfies 



From the work of [11] we know that the vectors U, V, W, X are propor- 
tional to a null tetrad that determines a self-dual metric. Indeed the tetrad 
is given by a a = / _1 V a , where V a = (U, V, W, X) for a = 0, 1, 2, 3 and 
f 2 = e (U, V, W, X), for e the four dimensional volume form dt Adx Ady Adz . 
In our case, f 2 = —h u and our line element is 

ds 2 = dt (hfy dy+h tz dz) + dx (h xy dy + h xz dz) + (h ty dy + h tz dz) 2 . (24) 

n>tt 
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3 The Formal Solution. 



We now construct, at least formally, the general solution to (|19|). Instead 
of working with this equation directly, it is helpful to define two functions 
A = h t , B = h x , and rewrite ( |T9"D in the equivalent form 



A t + A y B z - A z B y = 0, 



(25) 



A x = B t . (26) 
If we just viewing B as some arbitrary function, then the solution to (E5f) is 



A(t,x,y,z) = exp[/ dt 1 (B y (t 1 ,x,y, z) d z - B z (t 1: x, y, z) d y )\ ai(x, y, z), 



(27) 

where ai(x, y, z) is the value of A at t = as in (p0|). The exponential here 
is defined by its power series with the n'th term in this series being 

ft rti rtn-i 

/ dti / dt 2 . . . / dt n 

Jo Jo Jo 

[(fl v (ti)0* - B z {t x )d v ) ...{B y (t n )d z - B z (t n )d y )} ai {x,y,z). (28) 

We now must impose (p6j ) as a consistency condition on this solution. This 
gives us 



A(t,x,y,z 
exp 



t ptx 

dti / dt 2 (A xy (t 2 , x, y, z) d z - A xz (t 2 , x, y, z) d y ) 
o Jo 



ax(x,y,z). (29) 



Formally, this equation can now be solved iteratively. We can make successive 
approximations 

A® = ai{x,y,z), (30) 
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t 2 t 2 
A (1) = exp[((ta 0xy + —ai xy )d z - [ta 0xz + —ai xz )d y )}a l (x,y,z), (31) 

= exp [ J* dh jT* 1 dt 2 (A%j 9, - & ) ] ai(z, y, 2), (32) 

for n > 1. Then defining A = lim^oo gives the formal solution for 
v4Q. Integrating A with respect to t and imposing h(t = 0) = ao(x, y, z) then 



gives us a solution of ([191) 



Finally, we note that ( p5|) means that the quantity x, y, 2) is t inde- 
pendent, where y and 5 are defined implicitly by 

y{t) = y + f dhB^tux^ih),^)), (33) 
Jo 

z(t) = z - [ db x B v (t x ,x,y(t x ),z(t x )). (34) 
Jo 

This may be important if we were to look for action angle variables for the 
system. It also implies that A(t, x, y, z) = ai(x, y', z'), where the coordinates 
y', z' are defined by y(t, x, y', z') = y, z(t, x, y', z') = z. Thus the dynamics 
are characterised by a coordinate transformation in the y, z planeQ. 

4 Group Methods 

Several powerful techniques have been developed for the study of partial 
differential equations [12| . One of the most powerful is that of group analysis 
|j~3| , |i~4j1 . By studying the Lie algebra under which a given system of partial 
3 It is beyond the scope of this paper to show that the actually do converge to a 
well defined limit. 

An earlier version of this paper incorrectly stated that this transformation was area 
preserving. 
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differential equations is invariant, we can hopefully find new solutions to 
these equations. One method of doing so is to look for similarity solutions 
which are left invariant by the action of some sub-algebra of this symmetry 
algebra. This will reduce the number of independent variables present in 
the equation, possibly reducing a partial differential equation to an ordinary 
differential equation. However, such similarity solutions, by construction, 
will have some symmetries imposed upon them, so this method is not very 
useful if one is looking for the general solution to a system of equations. 

A more powerful method is to exponentiate the infinitesimal action of 
the Lie algebra into a group action, which takes one solution of the equation 
to another. However, even if this is possible, it is unlikely that the group 
action can be used to find the general solution to the equation from any given 
solution. 

Instead of attempting to find the symmetry algebra of (0) , it is easier to 
work with the equivalent system (|25|) and (|26|). We find that ( p5|) and ( p6|) 



admit a symmetry group defined by the infinitesimal generators 

£1 = f A d t - f x d B , (35) 

6 = (tg x + Bg A ) A d t + g A d x - g x d A - (tg x + Bg A ) x d B , (36) 

£ 3 = ktd t + kxd x + kydy, (37) 

&t = hd y - l v d z , (38) 

where / and g arbitrary functions of x and A, I is an arbitrary function of 
y and z, and k is an arbitrary constant. £3 just generates dilations, whereas 
£4 generates area preserving diffeomorphisms in the y — z plane. Although 
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^4 gives a representation of W m (modulo cocycle terms) [|T5| , they are really 
only coordinate transformations, so are not too interesting. However, we 
have two interesting symmetries, generated by £1 and £2- 

It is possible to exponentiate the action of £1 directly for an arbitrary 
function /. We find that if A (t, x, y, z) and B (t, x, y, z) are a solution of the 
system (|25|) and fl26|) then we can implicitly define a new solution, A and B, 
by 

A = A(t + f A (x, A),x,y,z),B = B(t + f A (z, A), x, y, z) + f x (x, A), 

(39) 

for any function / (x, A). Using this implicit form we can solve iteratively for 
the functions A and B given functions A, B and /. This means that given 



one solution of (19) , we can form an infinite dimensional family of solutions 
depending on that solution. For a given function g we can also exponentiate 
the action of £2, although its action cannot be exponentiated directly for a 
general function g. 

Although both ( |3T| ) and ( 153 ) give rise to infinite dimensional families of 
solutions from any given solution, they are not enough to derive a solution 
with arbitrary initial data from any given solution. 

If we compute the commutators of generators and ^(gj) for arbi- 

trary functions fa and gj we find that they obey the algebra 

[6(AUi(/ 2 )] = 0, (40) 

[£i(f)M9)] = Wa9* - f x 9A), (41) 

6(02)] = &(9lA92x - 9lx92A)- (42) 
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If we define a basis for transformations 

= Ca(x i+1 A m+1 ) (43) 
for a = 1, 2, where m and i are integers, then the above algebra becomes 

[(1)^™ (i) T n] = o, (44) 

[ (1) IT> (2) ^, n ] = ((™ + 1)0' + 1) - {n + m + l)) (1) ^+r' (45) 
[ (2) IT, (2) ^T1 = ((m + 1)0' + 1) - (n + 1)(< + l)) (2) ^r- ( 46 ) 
The algebra (fIBD is the algebra of locally area preserving diffeomorphisms 
which, modulo cocycle terms, is just the extended conformal algebra 
[ |T5|j . Thus (f44|) - ( f4"6|) represent some generalisation of W^. These are similar 
results to those found in |16| . 



We now note that equation (|T9|) can be derived from the Lagrangian 

f 1 1 

S = J d 4 x{-h 2 t + -h t (h y h xz - h z h xy )}. (47) 

The Hamiltonian is then 

H = \ L £x (7r ~ V h y hxz ~ h * h *v)) 2 - ( 48 ) 

where ti = h t + ^(h y h xz — h z h xy ) is the momentum canonically conjugate to 
h. We now define the Poisson Bracket of functionals of h and tt by 



The algebra ([H]) - fl46|) now reflects the fact that we have two infinite 
families of conserved quantities^] of the form 

I{f(x,A)) = ^ /(x,A)d 3 :c (50) 



3 We are assuming that we can ignore surface terms. 
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h(g(x,A)) = J^(tg x + Bg A )d 3 x. (51) 

These quantities all have vanishing Poisson brackets, i.e. they are in involu- 
tion. (The fact that these quantities are time independent comes from the 
conservation equations 

dtU) + fy(fB z ) - d z (fB y ) = 0, (52) 

and 

d t (tg x + Bg A ) + d x {g) - d y (tg xA B z + g A BB z ) 

+ d z (t g xA By + g A B B y ) = 0, (53) 

which follow from ( ^5|) and (|26|) .) 

5 Solutions 

We begin by looking for solutions that admit a triholomorphic Killing vector, 
£. This means the three complex structures, «T, are invariant under the action 
of £, i.e. J 1 = 0, where £ is the Lie derivative. Using the relationship 
between the complex structures and the vectors Vj given in Section 2 and 
the fact that £ is a Killing vector, we see that we require £gVj = 0. 

If d x is a triholomorphic Killing vector, this means that <9 X X = d x W = 
0, where X and W are as in ([16]) and (|T7|). This means h is of the form 
a(t,y,z) + xb(y,z) for some functions a and b. In terms of functions A = 
h t and B = h x this means that A = A(t,y,z), B = B(y,z), so ( p6| ) is 
automatically satisfied. If we take A(t = 0) = a\ (y,z) and B(t = 0) = 
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(y, z), it is straightforward to show that the solution to (B^) is then 

A(t,y,z) = exp{t((j) y d z - </> z d y )} a x (y, z), B (y, z) = (f)(y,z). (54) 
For given functions <p and ai it is straightforward to do the exponentiation, 



giving A explicitly. Using the exponentiated form of (|35| ) and (|36|), we could 
now use these solutions to generate new solutions which had some restricted 
x-dependent initial data as well. 

We can also consider metrics with a triholomorphic Killing vector, d z . 
This means we require 5 2 Vj = 0. In this case, we take h = —tz + g(t,x,y). 



We then recover the result [O] that g must satisfy the three dimensional 



Laplace equation g tt + g xy = 0. The general solution to this is known, and 
can be written in terms of two arbitrary functions ao(x,y) and ai(x,y). An 
almost identical reduction occurs if we take d y as a triholomorphic Killing 
vector. Again, using the symmetries (^) and fl3T)|), we can generate infinite 
dimensional families of new solutions, that in general have no Killing vectors. 

We note in passing that the solution corresponding to the multi-centre 
Eguchi-Hansen metric |T^| is 

A = - z + a £ sinh' 1 ( ( * ~ * j) _) , (55) 

*=i V2 v / (x-Xi)(y-^)/ 

B = _ a ' \fjt-U) 2 + 4(x-Xi)(y-rij ^ 
2 i=1 [x — Xi) 

where a is a constant. This is the only metric with a triholomorphic Killing 



vector that has a non-singular real (euclidean) section |L8|. 

At present, work is underway to complete a study of holomorphic Killing 
vectors and to see how the problem relates to the known results on such 
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metrics [O, [20[ . 



6 Conclusion 

We have shown how, at least formally, to construct the general complex 
metric with self-dual Riemann tensor. We have also studied the symmetry 
algebra of the system and found two infinite dimensional families of conserved 
quantities that have vanishing Poisson brackets. 

It should be emphasised that all the considerations here have been inher- 
ently local in nature, and we have imposed no sorts of boundary conditions 
on our solutions at infinity. If we were to look for metrics that are well defined 



globally, this would lead us to cohomological problems [21|, which appear to 
be best tackled using the twistor formalism [[|. 
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A The First Heavenly Equation 

Starting with (|25| ) and (|26|) , instead of looking on A as a function of £, x, y 
and 2 we take A as a coordinate and look on / = t and g = B as functions 
oi p = A, q = x, r = y, s = z. This transformation is well defined as long as 
A t 0. Inverting (|25|) and ( p6|) gives 

/? = -fe, (57) 

fr9s ~ fsQr = 1- (58) 

( |5"TD means we can introduce a function g, r, s) such that f = — Q p , g = 
Q q . (|58| ) then means that f2 must satisfy Q ps Q qr — Q pr Q qs = 1. Carrying 
out the same transformation on the line element (0), we find it becomes 
ds 2 = Q pr dp dr + Q ps dp ds + Q qr dq dr + Q qs dq ds. Thus we have recovered 
the Plebanski formalism. 
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